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Abstract. The concept of inner friction, by which a quantum heat engine is unable to follow 
adiabatically its strokes and thus dissipates useful energy, is illustrated in an exact physical 
model where the working substance consists of an ensemble of misaligned spins interacting 
with a magnetic field and performing the Otto cycle. The effect of this static disorder under a 
finite-time cycle gives a new perspective of the concept of inner friction under realistic settings. 
We investigate the efficiency and power of this engine and relate its performance to the amount 
of friction from misalignment and to the temperature difference between heat baths. Finally 
we propose an alternative experimental implementation of the cycle where the spin is encoded 
in the degree of polarization of photons. 
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1. Introduction 

The recent boosting interest in the study of the quantum counterpart of classical well-known 
heat engines such as Otto, Carnot, Stirling and Szilard ones 0], 2 3; 4] 5, 6, [7, 8j has 
been motivated both by the need of a fundamental understanding of the limits imposed by 
quantum mechanics on the thermodynamic performances of small devices (in terms of both 
efficiency and power output), and by the growing experimental ability to control various 
types of quantum systems with a high degree of accuracy. There have been, indeed, many 
proposals aimed at implementing thermodynamic transformations and cycles with many 
different quantum working substances, ranging from trapped ions to magnetic materials 0, 
with the prospect of building quantum heat engines, exploring the abilities and limitations of 
quantum machines in converting heat into work, and, on more general ground, to build a self 
contained description of thermodynamics in the quantum regime. As a specific example, an 
interesting proposal in this respect have been made for implementing a nanoheat engine with 
a single trapped ion, performing a quantum Otto cycle iflOll . Besides its specific applications, 
this is an important example as the quantum Otto cycle constitutes a useful test ground to 
study irreversibility in the quantum realm. 

Indeed, the cycle consists of two isochoric thermalization branches (with a fixed system 
Hamiltonian) and two isentropic branches in which the system is detached from the thermal 
baths and its evolution is generated by a parametric time-dependent Hamiltonian. Every 
practical realization of these latter transformations has to face the general problem of 
understanding and describing the (un-wanted) irreversible entropy production which can 
occur in non-ideal, finite time quantum parametric processes. This general problem has been 
variously analyzed through the use of fluctuation relations lfTTHT2l . and has attracted a lot of 
attention in recent years JT31 14J. 

In this paper, we explicitly address the study of the Otto cycle by focusing on the finite 
time case and discussing the implications of finite time transformations as opposed to ideal 
infinitely lasting ones. In this respect, in a series of papers, Kosloff and Feldmann [ 15] [T6l T71 
introduced the concept of intrinsic/inner friction, whereby the engine is never able to 
accomplish a frictionless adiabatic transformation and thus looses power. This concept has 
been then extended and applied to various contexts [ l2lfT8l 191. 

Inner friction is a fully quantum phenomenon, whose consequences are similar to those 
of the mechanical friction occurring when displacing a piston in compressing/expanding a 
gas in a classical thermodynamic setting. Its origin, however,, is completely different: when 
the external control Hamiltonian does not commute with the internal one, the states of the 
working fluid cannot follow the instantaneous energy levels, leading to additional energy 
stored in the working medium. Inner friction is thus associated to diabatic transitions, i.e. 
changes of populations which occur during the time dependent adiabatic (here referring to 
closed system) strokes if they are performed at finite speed. 

So far inner friction occurring in specific cycles and transformations has been analyzed 
by adopting phenomenological and physically motivated assumptions about the explicit time 
dependence of the population changes (e.g., in Ref. [15], a friction coefficient is introduced. 






Quantum Otto cycle with inner friction: finite time and disorder effects 


3 


giving rise to a constant dissipated power). Our treatment, instead, does not rely on any ad 
hoc assumption, but rather on the exact dynamics of the working substance. This is important 
because it has been shown, m, that inner friction is not only an indicator of irreversibility 
of a quantum process, but also a quantitative measure of its amount. It is therefore crucial to 
identify and highlight its role in the efficiency reduction of finite time cycles by analyzing the 
full quantum dynamics that produces it. 

In particular, we will explore the quantum friction arising from disorder within the 
sample playing the role of working substance. We will consider an ensemble of qubits in a 
setting in which their Hamiltonian parameters are not homogeneous, and connect the presence 
of these static errors to the appearance of friction and losses during the implementation of the 
Otto cycle. Explicitly, we provide a quantitative analysis of the amount of losses due to the 
inner friction as a function of the degree of disorder. 

Indeed, the performance of the heat machine are negatively affected by inner friction 
and cycle’s outputs such as extracted work, power and efficiency are gradually suppressed as 
disorder and friction increase. 


The remainder of the paper organized as follows. In Secs. |2.1| and |2.2| we introduce 
and review the concept of inner friction by focusing on the particular case of a spin system 
in presence of misalignments and disorder, which will then be of interest for the rest of 
the present paper. In Sec. |2.3| we introduce the quantum Otto cycle and its constituent 
transformations specifying the assumptions about the model we use to describe the working 
substance. In Sec. [3] we present and discuss our main results, while in Sec. [4] we propose a 
feasible experimental implementation of the quantum Otto cycle in order to test our findings. 
Finally, Sec. [5] is devoted to some concluding remarks and to a discussion of possible future 
developments. 


2. Model and methodology 

In this section we introduce the model, and give a possible explanation of the origin of inner 
friction. We then introduce the quantum Otto cycle (QOC) and the figures of merit through 
which the cycle will be characterized. 

2.1. Misalignment and disordered samples 

In order to understand what we mean by losses and friction in a closed quantum system, and 
in particular in the case of one qubit, let us focus on the dynamics generated by a Hamiltonian 
of the form: 

H(X(t)) = ~^ a z + A(f)(cos 9a z + sin 9cr x ). (1) 

The analysis reported here applies to the general case of qubit dynamics ([!]) and in the 
following we will consider the case of a spin interacting with a magnetic field (Sections 2 
and 3) as well as of a qubit encoded in photon polarization (optical implementation in Section 

4). 
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Figure 1 . Misalignment effects of a sample of spins interacting with an external field due to 
(a) disorder within the sample or (b) lack of uniformity of the magnetic field (e.g. because of 
the finite length of the coil). 


An adiabatic transformation is obtained by the unitary time evolution generated by the 
Hamiltonian (jT|), with a linear driving of the external field at a fixed rate X(t) = au> 0 t/2, which 
we allow to be misaligned by an angle 6 with respect to the static field ui 0 . The misalignment 
affects the energy spacing as well as the eigenstates and the populations. 

We assume that at t = 0 the qubit in a thermal state at inverse temperature 3. For a very 
slow driving, ideally taking an infinite time to complete the transformation in the quantum 
adiabatic regime, the qubit populations would remain unchanged while the the energy spacing 
increases/decreases and the system remains in a thermal state with a lower/higher temperature. 
The same occurs in absence of misalignment, 9 = 0 in Eq. W. as in this simple scenario the 
adiabatic transformation 

H = l o(t)a z , with u(t) = — + A (t) (2) 

reduces to a compression/expansion of the energy spacing of the qubit, thus preserving the 
initial thermal populations even in presence of fast driving. 

Interesting dynamical and thermodynamical implications arise, instead, when consider¬ 
ing deviations from the limit of perfect alignment (6 f 0). This is the case that we are going 
to explore in this work in order to characterize inner friction and its effects on the efficiency 
of quantum thermal machines, resulting from the simultaneous presence of the static field oj 0 
and of the misaligned time depending part A (£). Our aim it to apply this analysis to an en¬ 
semble of (independent) spins, considering some degree of disorder and looking at average 
effects across the sample. In particular, this can correspond to different situations as repre¬ 
sented in Figure [TJ A condensed system on a lattice, with embedded magnetic dipoles having 
disordered orientations, can be modelled by randomly oriented spins with tilting angles 9i 
(i = 1,2,...) with respect to the direction of a uniform external field. We assume that the 
distribution of the spin orientations in the sample is given by a function G(6). Alternatively, 
all sample dipoles could be perfectly aligned in an ordered configuration but the inner friction 
could be due to inhomogeneity of the external fields in space (Figure [!})). The field orientation 
across the sample would be given, in this case, again by the function G(6). 

2.2. Inner friction and irreversibility 

In order to have a simple physical picture for the behavior of our quantum machine, let us 
first consider the simpler case of a driven quantum two level system undergoing the unitary 












Quantum Otto cycle with inner friction: finite time and disorder effects 


5 


dynamics generated by a parametric time dependent Hamiltonian H[X(t)]. If the parameter 
A (t) changes slowly enough (in the sense of the quantum adiabatic theorem, [20, 21]) the 
system evolves without its energy population ever changing at all, even if the instantaneous 
energy eigenvalues and eigenstates do change in time. If the system has been prepared in 
equilibrium with a thermal bath, which is then removed, such an ideal adiabatic parameter 
change keeps the system in an equilibrium state at every stage. In particular, if the parameter 
A gets back to its initial value after some time, the final result is that the system is brought back 
to its initial state. On the other hand, if the cycle is performed in finite time, the final state of 
the system will differ from the equilibrium state it started off because non-adiabatic transition 
have taken place Il22ll . The difference between the two states, if properly quantified, can be 
regarded as a measure of the deviation from an ideal adiabatic transformation. The quantum 
non-adiabaticity has the same effects as friction has in a classical context: an extra energy is 
needed to complete the process (indeed, the work done in the ideal adiabatic is always smaller 
than the one performed in finite time, see lf2Tl [T2l0 - which is then dissipated if the system 
equilibrates at the end of the process. 

With this picture in mind, let us now address the dynamics generated by Hamiltonian ([!]) 
on an initial thermal state given by p 0 = exp{— fiH(X(0))}/ Tr[exp{—/5i7(A(0))}] where ft 
is the inverse temperature in units of the Boltzmann constant. By changing X(t) very slowly 
from A(0) at t = 0 up to X(tf) = X* at t = tf and then going back from A* to A(0) the system 
will be brought back into its initial state. To discuss what happens in the general case, namely 
when these changes are performed at finite rates, we consider the following protocols: 


p2 < - Pi 

Ub{PJb) 

The forward protocol, defined by the unitary operator Up ft), t F ) = Hix i^ T )) dT 

(T being the time ordering operator), is generated by the Hamiltonian in Eq. ([!]) such that 
A pit) = O'pCcVji/2. It takes the initial density matrix p 0 to p\ = Up ft).tp) p 0 Up (0,t F ). The 
backward protocol U B (0, tpft) = Te~ l l° B H{\ B {j))dT a g a j n generated by the Hamiltonian in 
Eq. 0 where now X B {t) = X F {t F ) — a B co 0 t/2 with the condition that X B (t B ) = X F (0). 
This consists just in ramping up and down the field X(t) with different rates a F and a B , 
respectively. 

In order to characterize the above protocol, we first look at the time dependent 
polarization defined as n(t ) = Tr[p(t)H(t)]/u(t) where u(t) is the energy level spacing 
at time t for both the forward and backward protocols. The result is shown in Fig. |2(a) , 
where we notice that finite-time evolution introduces deviations with respect to the quantum 
adiabatic case as expected. Moreover, as it can be seen in Fig. |2(b)[ by applying the forward 
and backward protocols defined above, the system does not get back to its initial state, but 
reaches a different polarization (green line) at the end of the protocol. This already gives 
a quantitative indication that finite time control leads to an irreversible behavior. Here, we 
use the word “irreversibility” in the thermodynamic sense: because of the occurrence of 
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Figure 2. The left panel shows the time dependence of the polarization n(t) for a finite-time 
adiabatic transformation (solid blue line) compared with the case of an ideal quantum adiabatic 
one (dashed red line). We used ap = 10 _4 wo- In the right panel, we compare the forward 
evolution (lower time axis from left to right) with the backward one (upper time axis from 
right to left) by displaying the time dependent polarization in both processes. The parameters 
used are aptp = cub^s = 15 and ap = cub = ^o- In both figures 9 = n/5. 


non-adiabatic transitions, the system is driven out of the manifold of equilibrium states and 
application of the same protocol in reverse does not bring it back to the initial state. 

A more precise way of quantifying the irreversibility of such a transformation is through 
the distance of the final state from the initial one, expressed in terms of the relative entropy 
D(p 2 \\p 0 ), where p 2 = Up(0, t D ). As shown in lfl2l . this quantity has a well 

defined thermodynamical interpretation as it precisely gives the non-adiabatic part of the work 
performed on the system by the driving agent, i.e. the inner friction. 

Indeed, for an adiabatic transformation, the quantum relative entropy between the actual 
final state and the ideal thermal equilibrium one is proportional to the difference between 
the work done on the system during the parametric change and the same quantity taken in 
the infinitely slow limit. This is precisely the definition of the inner friction, hereafter called 
Wfric, m. Furthermore, the same quantity is linked to the generation of extra heat; that is to 
say, to the irreversible production of ‘waste energy’. This extra energy is exactly the energy 
that needs to be dissipated if, at the end of the protocol, we were to thermalize the system to 
the initial temperature. Specifically, the following relations hold: 

— (3Q(p2 Po) = fiWfric = D(p 2 \\po) ( 4 ) 

where Q ( p 2 —>■ pf) is the heat the system takes to thermalize at the initial inverse temperature 
f3. This is what we shall refer to as inner friction in the following. 

The inner friction for the time evolution described above is reported in Fig. [3j where we 
can clearly see that when both transformations are either very slow (quantum adiabatic case) 
or very fast (‘diabatic’ or sudden case), at the end of the protocol the system is found to be in 
(or very close to) its initial state. For finite time transformations, however, the system does not 
get back to its initial state. From a dynamical point of view, this is not surprising; however, if 
interpreted from a thermodynamical perspective, this fact suggests that transformations done 
in finite time are, in general, irreversible ones. 
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Figure 3. The quantum relative entropy between the state at the end of the backward step and 
the initial thermal state po. Notice that inner friction is very close to zero both for ap,B —> 0 
and op B —> oo. The inverse temperature characterizing the initial state is taken to be equal to 
the energy spacing cvq. Here 9 = 7r/5 and a pip = a pip = 15. 


2.3. Model for the quantum Otto cycle (time scales assumptions) 

The Otto cycle is the simplest cycle for our purposes as it allows for a clear separation 
between dissipative steps (thermalization processes, in contact with a thermal bath) from non- 
dissipative ones (in which work is done or extracted), as opposed, for instance, to the Carnot 
cycle, which contains two isotherms in which one has to perform (extract) work while the 
system is attached to a thermal bath. This separation will be very useful in order to identify 
finite time effects on the single adiabatic transformations and thus on the total cycle. 

The quantum version of the Otto cycle is the composition of two adiabatic 
transformations, in which the systems evolves unitarily, and two isochoric branches 
corresponding to thermalization in contact with a hot (and, respectively, a cold) heat bath 
at temperature 8f l (j3f l ). 

In the next subsections we better specify the assumptions we employ to describe the 
different branches of the quantum Otto cycle and the physical quantities we investigate to 
characterize it. The ideal Otto cycle is represented by the dashed (yellow) rectangle in Fig. |4} 
The blue line, instead, describes a finite time cycle, in which the end points of the adiabatic 
strokes are moved towards larger values of n (which just means that there is more population 
than expected in the excited states) because of the presence of inner friction. 

2.3.1. Adiabatic transformation - As already mentioned above, the adiabatic transforma¬ 
tions can be described by the unitary operator generated by the Hamiltonian in Eq. O- For 
simplicity, in the following we will consider the case where the two adiabatic branches (1 —>■ 2 
and 3 —>■ 4) last equally long, namely tp = tp = t 0£ j and have the same rate of change for the 
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Figure 4. Representation of the Otto cycle in a parameter space in which the horizontal axis 
gives the instantaneous energy spacing between the eigenstates of the Hamiltonian ([]]), while 
the vertical axis gives the polarization. The solid blue line is an example of finite-time Otto 
cycle with parameters: 6 = at = uo^t = 0.5513, fa = uo^ 1 and fa = fa/2. The 
dashed yellow line corresponds to an ideal (infinite time) Otto cycle. The two red lines are 
the isotherms in this plane. They include the two (very fast) branches in which the system 
equilibrates in contact with baths at inverse temperatures fa and fa, respectively. 


field Oip — OLft. 


2.3.2. Isochoric transformations - For the isochoric transformations, we assume perfect 
thermalization at the given temperatures 3f l , (hot and cold, respectively). To study the 
relation between inner friction and the finite time of the adiabatic branches, we will assume 
that perfect thermalization is achieved very fast with respect to all other time scales, and also 
that the isochoric branch will be assigned with a fixed short time duration (to be eventually 
neglected) with respect to adiabatic ones but long with respect to thermalization time of the 
system: 

1~therm < ^- 7'iso 1~ad (5) 

where r t h erm , r lso and t 0( i are the typical time scales for the thermalization process, isochoric 
and adiabatic transformations, respectively. 

2.4. Figures of merit 

In order to characterize the quantum Otto cycle, we will look at the extractable work, W ex , at 
the power V, at the efficiency rj, and at their averages over disorder. 
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To properly define these quantities, let us start by defining the work done on an adiabatic 
branch as: 

W = Tr [H f p f \ - Tr [H iPi \ (6) 

where Hi (Hf) and pi (pf) are the Hamiltonian and the density matrix of the system at the 
beginning (end) of each transformation. In particular, both adiabatic transformations start 
with a Gibbs-like state since we assume perfect thermalization to occur at the end of each 
isochores. In the adiabatic transformations, the work defined in Eq. ([6]) does coincide with 
the first moment of the work distribution for closed but non-autonomous systems Il25ll . On 
the other hand such a work distribution allows to define a fluctuation relation and thus its 
moments have a clear thermodynamical meaning. 

In the isochoric branches, we have that the initial and final Hamiltonians are the same 
and the final state pf is thermal, and thus diagonal in energy eigenbasis. The amount of energy 
exchanged between the reservoir and the system in each isochoric transformation is given by: 

Qiso = Cl (p[ f) - Pi°) + e o ( Po f) - Po ’) = W (po } - Po /} ) • (7) 

Thus, the energy absorbed from the bath equals the energy spacing oj = e\ — e 0 , times 
the change in the population of the lowest energy state (we denoted the ground and excited 
state populations as p 0 andpi, respectively). 

Since the change of the total internal energy along the cycle vanishes, the total work 
done on the system is given by W tot = — ( Q h + Q c ) where (f, ((f ) is the amount of energy 
exchanged with the reservoir at inverse temperature ff L (fie), given by Eq. (|7|) for the isochores 
2 —> 3 (4 —> 1). The first quantity we will use to characterize the cycle is the extractable work 
W ex = -W tot = (Qh + Qc) given by the relation: 

w ex = (u2(Po ] - P ( o 3) ) + Ui(p ( 0 4> - P ( 0 1] )) , (8) 

where uj,. = ) is the energy level spacing of the Hamiltonian at point k = 1,2 in 

the cu-n diagram of Fig. |4} 

For the quantum Otto cycle, and by means of the definition of n, we can then write the 
following condition 

- n (4) ) < w 2 (n (2 ) - n (3 )), (9) 

ensuring that the work extracted is strictly positive and we are actually using the engine to 
perform work. This is in agreement with Carnot theorem, as shown in Ref. Ifl5ll . and for ideal 
quantum adiabatic branches (i.e. infinitely slow transformations) it reduces to IfJ f. > coi/u >2 
where c<Ji j2 arc the energy spacings at the end of quantum adiabatic branches [23J. Notice 
that, in this quantum adiabatic limit, a breakdown of the positive work condition coincides 
with the saturation of Carnot inequality. 

Let us now define the power and the efficiency of a cycle on a single qubit as: 
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V{Tad,0,P h /P c ) 

r](T ad ,6,f3 h /f3 c ) 


W4 


W4 


tp + ts + T~iso ‘^■T~ad 

ex Q c 

Qh Qh 


Tit 


( 10 ) 

( 11 ) 


where we have made explicit the dependence upon the angle 9 in both V and rj and we 
explicitly included r iso in the above definitions, despite the assumption that it is small 
compared to the characteristic times scale of the cycle, just to avoid that V is ill-defined 
in the formal limit T ad —» 0. 

In addressing the disordered case we average all quantities over 9 assuming it to have a 
Gaussian distribution G a (9) with zero mean and variance a 2 . The averaged extractable work, 
power and efficiency are thus given by: 


W ex (T ad ,/3 h //3 c ,o) 


V{r ad , fi h /fi c ,<r) = 
rj(r ad ,p h /l3 c ,o) = 


= / G a (9)W ex (r ad ,9ffi h /ffi)d9 

J 0 

r G,(9)V(r adM) dO 


G a (6)ri(T ad ,9,f3 h /f3 c )d9. 


3. Results and discussion 


( 12 ) 

(13) 

(14) 


In this section we characterize the QOC by looking at the extractable work W ex , its power V 
and its efficiency r/, paying particular attention to the role of the inner friction in limiting the 
performances of such an heat engine. In the first subsection we will look at these figures of 
merit for different sets of parameters of our model-system. In the second one, we study the 
behavior of the efficiency at maximum power, t](Vmax)- 


3.1. Extractable work, power and efficiency 


In what follows we will study the performance of the QOC in various cases of equal rates for 
both of the adiabatic branches. In Figure |5(ajl|5(c)| and |5(e)| we plot W ex , V and ?/ as function 
of the total time of the cycle t fot for different values of the misalignment angle 9 at a fixed 
value of the ratio between the temperatures of the hot and cold reservoirs ffi /= 0.5. It can 
be seen that the extractable work becomes negative if the t tot exceeds a maximum time 
which is a function of the misalignment 9. This means that if the cycle lasts too long we are 
actually doing work on the system. Moreover there exists a critical value of 9 such that the 
extractable work is negative for any value of t tot . Under these conditions, the cycle is not a 
heat engine but rather a refrigerator, which uses external work to cool the cold reservoir. An 
analogous behavior is shown by power and efficiency. 

In Figure |5(b)[ |5(d)| and |5(f)[ we show the dependence of W ex , V and r\ on the total 
time of the cycle, for a fixed misalignment angle 9 = 7t/ 5 and for different values of the 
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(a) Extractable work ( W ex ) vs. attot for (b) Extractable work (W ex ) vs. at to t for various 
different values of 0 fa/fa 



OttroT 

(c) Power (V) vs. at tot for various of 0 



O'tror 

(e) Efficiency ( 77 ) vs. at tot for various 9 




(f) Efficiency ( 77 ) vs. at tot for various fa/ fa 


Figure 5. Extractable work, power and efficiency as functions of at to t • In Figure |5(a)| |5(c)| 
and 5(e) we fix the temperature ratio fa/fa = 0.5 and vary the misalignment 0 , starting from 
0 = 0 for the highest (red) plot, then considering 0 = 7 r/ 5 , 7 r /2 and finally 0 = tt for the 
lowest (yellow) curve. On the other side, in 5(b) 5(d)| and [5(f)| we fix the misalignment to 
0 = 7r/5 and vary fa/ fa, which again going from the top-most red plot down to the yellow 
curve takes the values fa/fa = 0.01,0.31,0.51,0.71. 


ratio fa/ fa. We can see that as the ratio increases the extractable work increases too. This 
is something which is expected, nevertheless we can clearly see that the finiteness in time 
of the cycle introduces again negative works for t tot > tu(fih/ftc)- This is again due to the 
generation of inner friction, which comes along with the finite time condition. 

Inner friction is explicitly shown in Figure [6} where the sum of the friction produced 
in the two adiabatic strokes is shown as a function of the total cycle time for various 
misalignment angles 9. As discussed above, the case 9 = 0 is very special as no friction 
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Figure 6. Inner friction accumulated in the cycle as a function of the total operation time for 
different misalignments 0, at /3h//3 c = 0.5. 


is generated, whatever rate of variation is considered for the driving field \(t). On the other 
hand, Wf ric increases with increasing the tilting angle 0 and decreases with decreasing the 
driving rate a. The behavior of Wf r i C should be compared to that of W ex shown in Figure [5j 
the more friction is present, the less work can be extracted from the engine. 

So far we have considered QOCs with given values of the misalignment angle 6. We 
now consider the effect of disorder and assume that 9 is a Gaussian random variable with mean 
value 9 = 0 and variance a 2 . In Figure 7(a), 7(c) and 7(e) we show the behavior of extractable 


work, power and efficiency for different values of the variance and given temperature ratio 
(fih/fic = 0.5). We can see that at a given total time t tot , the best performance is always 
obtained with sharper distributions (smaller a). Thus, if the disorder of the system grows, the 
capability of the latter of providing work and of doing it in a more efficient way decreases. 
Again we mention the fact that there exist a maximum total time above which the QOC 
is not an heat engine anymore. We also notice that even a small disorder has quite dramatic 
effect in reducing the efficiency for long enough times (upper curve in Fig |7(e)[ ). 

In Figure |7 (b)[ |7(d)| and |7 (f)1 we plotted the behavior of the same quantities for different 
values of the ratio fih/fic at a given variance a 2 = 0.1. Again, all of the quantities increase as 
the difference in temperatures increases. 

We plotted all of the quantities as a function of the total cycle time t tot because, from 
an operational point of view, this is the quantity one can control once the working substance 
is prepared and the stage is set for the thermal machine to operate. However, it has to be 
mentioned that some care should be payed when comparing the values of the efficiency at 
different operating times. Indeed, since uj varies with time, each different t fot gives rise to a 
different value of the final frequency (called ui 2 , above) at which the isochoric 2 —» 3 takes 
place, see Figure (|4]). The ideal cycle with infinitely slow adiabatic branches, corresponding 
for us also to the absence of misalignment (9 = 0) and shown for comparison in each plot in 
Figure (|5j) and (|6j), has efficiency riueai = 1 — The dependence of r/ideai on implies 
that the efficiency r\ should be compared with a different r/id ea i at each different t tot . To avoid 
any confusion in this respect, and to better display the role of finite-time induced friction in 
the machine performance, we show this comparison in Figure ([8]), where the ratio r//r/ M | ea ; is 
displayed as a function of the operating time. Once the efficiency is renormalized in this way, 
its residual dependence on t tot can be fully ascribed to the presence of inner friction. 
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(a) Average extractable work (W ex ) vs. attot (b) Average extractable work (Wex) vs. at to t 
for various cr for various Ph / p c 



0 1 2 0 1 2 

OttroT Cx£ T ot 


(c) Averaged power (V) vs. at tot for various cr (d) Averaged power (V) vs. at tot for various 

Ph/Pc 
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(e) Averaged efficiency (jf) vs. at to t for various (f) Averaged efficiency (rj) vs. at to t for various 

O' Ph/Pc 


Figure 7. Averaged extractable work, power and efficiency as functions of attot • In Figure 
7(a) 7(c) and 7(e) we consider a fixed value for the temperature ( Ph/Pc = 0.5 ) and vary a. 


Different colors refer to different gaussian widths: the red plots correspond to a 2 = 0.01, the 
orange ones to cr 2 = 0.5, light orange ones to a 2 = 1 and finally the (lowest) yellow curves 
refer to a flat distribution. In 7(b) 7(d) and 7(f) we take cr 2 = 0.1 and vary Ph/Pc> which, 
going from the red plots to the yellow ones, takes the values Ph/Pc — 0.01,0.31,0.51,0.71. 
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atjoT 



Figure 8. Left: renormalized efficiency T]/r]id ea i as a function of the total operation time for 
different misalignments 9 = 0, n/10, n/5,2ir/5, at dh/ dc = 0.5. Right: averaged efficiency 
fj normalized with respect to the ideal efficiency obtained at 9 = 0. The average is taken 
over gaussian distributions with variances, a 2 = 0.1, a 2 = 1, and over a flat distribution, 
respectively. For all of the plots, we fixed the temperature ratio dh/dc = 0.5. 



Figure 9. Relation between averaged power and efficiency, V(t) and 77 (f) at the same time 
parameter attot • In 9(a) we choose dh/dc = 0.5 and vary < 7 , which, from the outer to 
the inner curve takes the values a 2 = 0.01,0.05,0.1,0.5,1,10. In 9(b) we fix a 2 = 0.1 


and starting from the outer to the inner curve we consider the increasing temperature ratios 
dh/dc = 0.02,0.03,0.04,0.05,0.06,0.07. 


3.2. Efficiency at maximum power 


Let us now consider the relation between rj and V (see Figure[9]) and then extract the value 
of the efficiency at maximum power t](Vmax)- Two sets of these data are reported in Tab. 
3.2|and |L2l They refer to averaged power and efficiency considering in the first case various 
temperature ratios at a fixed width, and the other way round for the second case. The effects 
of disorder and temperature difference continue to stand: our analysis provides larger values 
of (averaged) power and (averaged) efficiency at maximum power for smaller and smaller 
<7 and for larger and larger dh/ fie with the following one-sentence summary: We obtain 
considerably larger values of t](Vmax ) when the temperatures ratio is large and for very 
picked misalignment G distributions, that is, when the inner friction is smaller. 
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a 2 

at MAX 

V max/ o? 

max) 

0.01 

0.0882 

0.0439 

0.0775 

0.05 

0.0882 

0.0429 

0.0758 

0.1 

0.0882 

0.0418 

0.0737 

0.5 

0.0771 

0.0334 

0.0519 

1 

0.0340 

0.0253 

0.0340 

10 

0.0.220 

0.0027 

0.0220 


Table 1. Efficiency at maximum power at /3^//3 c = 0.5 and for different values of the Gaussian 
bell’s width a. The optimal total cycle time, atf/f x , is the one for which V attains its 
maximum. 


Ph 

a t MAX 

V max/ o? 

vi/P max) 

2.1 

0.175 

0.0761 

0.1420 

3.1 

0.125 

0.0635 

0.1056 

4.1 

0.1 

0.0517 

0.0862 

5.1 

0.075 

0.0406 

0.0660 

7.1 

0.05 

0.0208 

0.0447 

9.1 

0.025 

0.0045 

0.0224 


Table 2. Efficiency at maximum power for a 2 = 0.1 for different temperature ratios Ph/ Pc- 
The maximum power V max = is obtained for the times reported in the 

second column. 

4. Experimental implementation with an optical set-up 

In this Section we propose one possible experimental set-up by means of which it is possible 
to realize the Otto cycle discussed so far and test our findings. The Otto cycle is made up 
of two different types of branches, namely adiabatic and isochoric transformation. Thus the 
proposed set-up has to be able to implement both of them. 

The physical system we have in mind is an optical one and in particular we propose to 
encode the qubit into the polarization degree of freedom of a single photon. In the following 
we address the implementation of the two types of branches separately, stressing the key 
points for both of them. 
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4.1. Implementation of the adiabatic transformation 

The adiabatic branch in the Otto cycle is achieved by means of a time dependent effective 
magnetic field and its time evolution is given by the unitary operator: 

U(t) =Te~ l ti dT S{ - T) ' 3 (15) 

For a fixed t = t* the above operator can be written as a rotation in the Hilbert space of 
the qubit using the Euler decomposition as: 

(7(f) = (16) 

This expression is helpful for our purposes because the single rotations appearing in it can be 
easily implemented in an optical setup as rotation of the polarization degrees of freedom of a 
single photon. 

Therefore by encoding the qubit into the polarization degree of freedom of a photon and 
in particular by choosing the basis {\H) , |V)} of horizontal and vertical polarization, we can 
perform the wanted rotations by means of properly chosen phase retarders. 

4.2. Implementation of the isochoric transformation 

The isochoric transformation requires more care. By definition it amounts to attaching the 
system to a thermal bath, which makes the system thermalize into a Gibbs-like state. The 
latter is characterized by a density matrix p with no coherences between different eigenstates 
of the Hamiltonian of the system whereas the diagonal ones are given by the Boltzmann 
factors e -/3en /(e -/3e ° + e _/3ei ), where ft = 1/T C or ft = 1/Th is the inverse temperature we 
want the state to thermalize at and e n are the eigenenergies of the final Hamiltonian of the 
adiabatic transformation preceding the isochoric one we are addressing. In the case of a qubit 
the thermalization process leads to a final state which is diagonal in a given basis (determined 
by the form of the bath-spin coupling) and whose populations are related to the temperature 
of the thermal bath by the relation 

ft = T~ 1 = -log^- -<p { 0 f) < 1, (17) 

where, p^ is the population of the lowest (n = 0) state of the qubit after thermalization has 
occurred. 

In order to implement an isochoric transformation we propose to exploit the experimental 
set-up used in Ref. [24]. The idea is to exploit the spatial degrees of freedom of the photon as 
an effective bath for its polarization. The coupling between the two is achieved by exploiting 
the birefringent property of a quartz plate. The effect of the latter on a photon passing through 
it is to phase-shift the horizontal and vertical component of the polarization by an amount 
proportional to the number of photons per mode. Once the spatial part of the photon is traced 
out, the dynamics of the polarization turns out to be driven by the following dynamical map 
between an initial state pi to a final state pf, which describes decoherence: 

Pi ^ Pf = 2 ((1 + z ) Pi + (1 — z ) ®z Pi ftz) ) 


( 18 ) 
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where the parameter z can be tuned from z — l (identity map) to z = —1 (complete 
decoherence, namely the final density matrix has vanishing off-diagonal terms). Because 
of our assumption of complete thermalization we will always assume z = — 1. 

We can thus exploit this mechanism in order to engineer thermalization in the following 
way. Let us assume that the inverse temperature of the bath we want to mimic is f3. 
Through Eq. ( fT7| ) we can determine the population of the lowest energy level after the system 
completely thermalized. Let us write the initial state (which in turn correspond to the final 
state of the adiabatic transformation preceding the isochoric one) as: 


Pi 




O z T h 'x (7 X T by (7 y 


(19) 


Since the decoherence mapping in Eq. ( [18] ) has the effect of making the off diagonal elements 
vanish we first need to perform a rotation on the initial state pi to turn it into a state of the 
form: 


p'f=o 1+ o -p 0 


If) 


+ b' x e x + by o y 


( 20 ) 


where p ( f > is calculated through relation in Eq. (|17j). The application Eq. (1181) has now the 
effect of making b' x = b[ ' = 0, thus leaving us with the desired state: 


Pf = \ 1 + Q -Po /} ) Oz ( 21 ) 

It is easy to see that in order to get from the state in Eq. ( fl9| ) to the one in Eq. ( |20| ) we can apply 
a specific rotation, which has to be chosen by taking into account both states. For instance in 
the case b y = 0 and p\P < p$ <1/2, that is, if we are ‘heating’ our system, the right rotation 
to perform is: 


R 


X 




= e 


Ox 

2 




( 22 ) 


with cos(d,.) = (l — 2 p ( fj (l — 2 pI^J . For b y / Owe have solutions for cos (0 y ) only if 

/ 2 \ 1/2 

-l/2<l/2-p ( P <-fl/2-p^) +bl J . 


5 . Conclusions 

We have worked out an exact dynamical model with which we discussed the performance of 
a Quantum Otto cycle in presence of inner friction. With respect to previous approaches to 
the same problem, we obtained the growth of polarization in the adiabatic branches of the 
cycle without any ad-hoc assumption, but rather by following the dynamics generated by the 
system Hamiltonian. In this way we have been able to deal with the irreversibility of such 
transformations when they are performed in a finite time and so to better characterize the 
whole cycle. Finite time evolution leads to a decrease of values of thermodynamical figures 
of merit for the heat engine, and we have concentrated on the extractable work, the power and 
efficiency which are partially quenched if inner friction is present. 
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The friction is related to the non-commutativity of the system and control Hamiltonian, 
due to some misalignment between the internal and control magnetic field axes. After 
explicitly studying its effects for a fixed (and controlled) case, we turned to the more realistic 
case in which such a misalignment is an unwanted side-effect of the lack of control in the 
system, ultimately due to the presence of disorder in the sample or inhomogeneity of the 
magnetic field. The limiting cases are the completely ordered and disordered samples; in the 
first one, we obtain an ideal quantum Otto engine with efficiency given by r] u i e ai = uj 2 ; 

while for the completely disordered case we showed that no positive work can be extracted 
from the system, which cannot behave as an heat engine at all. In the more general case of a 
finite-width Gaussian distribution of tilting angles 9, describing a given degree of disorder as 
quantified by its variance, we obtained a quantitative description of the efficiency reduction 
due to the disorder-induced inner friction. 

Finally, we proposed an optical experimental implementation of such Otto cycle using 
the effective polarization qubit of a photon, whose thermalization can be obtained by the 
coupling with the spatial degree of freedom in a birefringent crystal. 
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